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Abstract 

The divergent large-order behaviour of the perturbative series relevant for 
the determination of a s from r decay is controlled by the leading ultraviolet 
(UV) renormalon. Even in the absence of the first infrared (IR) renormalon, an 
ambiguity of order A 2 /m 2 is introduced. We make a quantitative study of the 
practical implications of this ambiguity. We discuss the magnitude of UV renor- 
malon corrections obtained in the large- iVy limit, which, although unrealistic, is 
nevertheless interesting to some extent. We then study a number of improved 
approximants for the perturbative series, based on a change of variable in the 
Borel representation, such as to displace the leading UV renormalon singularity 
at a larger distance from the origin than the first IR renormalon. The spread of 
the resulting values of a s (m 2 ) obtained by different approximants, at different 
renormalization scales, is exhibited as a measure of the underlying ambiguities. 
Finally, on the basis of mathematical models, we discuss the prospects of an 
actual improvement, given the signs and magnitudes of the computed coeffi- 
cients, the size of a s (m 2 ) and what is known of the asymptotic properties of 
the series. Our conclusion is that a realistic estimate of the theoretical error 
cannot go below 5a s (m 2 ) ~ ±0.060, or 8a s (m 2 z ) ~ ±0.006. 
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1. Introduction 



The possibility of measuring a s from T-decay has been extensively studied in a 
series of interesting papers, in particular by Braaten, Narison and Pich The 
relevant quantity is R T = T(r — ► v T + hadrons)/r(r — ► v T + lu), with I — e, At 
present || the ALEPH collaboration finds 

( R r) eX p = 3 ' 645 ± °- 024 - (1-1) 

Based on this result, it is argued, if the formalism of QCD sum rules is assumed, 
according to SVZ ||, that 

a s {m 2 T ) = 0.355 ±0.021 (1.2) 

and finally 

a s (m 2 z ) = 0.121 ± 0.0016(exp) ± 0.0018(U) = 0.121 ± 0.0024. (1.3) 



Given that m T is so small, this determination of a s (m|) appears pj-|T0"|] a bit too 
precise! 



In defence of this method |TT[ one can certainly point out that R T has several 



combined advantages. Dropping some inessential complications, R T is an integral in 
s of a spectral function R(s) which is the analogue of R e + e -(s) but for the case of 
charged weak currents. Thus, first, R T is even more inclusive than R e + e -(s) and one 
expects that the asymptotic regime is more precocious for more inclusive quantities. 
Second, one can use analyticity in order to transform the relevant integral into an 
integral over the circle \s\ = m 2 This not only gives some confidence that the 
appropriate scale of energy for the evaluation of R T is of order m T , but also shows 
that the integration over the low-energy domain helps very much in smearing out 
the complicated behaviour in the resonance region. Also important is the presence 
of a phase-space factor that kills the sensitivity of the spectral function near Re s = 
m 2 ., where there is a gap of validity of the asymptotic approximations due to the 
vicinity of the cut singularities and also to the nearby charm threshold. On the circle 
\s\ = m 2 , asymptotic formulae should be approximately valid for the correlator. The 
perturbative component of R(s) is known up to terms of order a s (m 2 ) 3 [|12| . One can 
hope to get some control of the non-perturbative corrections by using the operator 
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product expansion and some estimate (either experimental or by some model) of the 
dominant condensates, in the spirit of the QCD sum rules ||. 

This series of virtues of R T is indeed real but would not be sufficient in itself to 
justify the precision on a s (m 2 ) which is claimed. The real point is that no corrections 
of order 1/m 2 are assumed to exist. The fact that there is no operator with the cor- 
responding dimension in the short distance expansion is not sufficient, because there 
could be non-leading corrections in the coefficient function of the leading operator. 
We think it is a fair statement that there is no theorem that guarantees the absence of 
A 2 /m 2 terms in R T in the massless limit; no theorem that proves that terms of order 
A 2 /m 2 cannot arise from the mechanism that generates confinement. But even if in 
principle the above theorem would exist, still, in practice, there would be ambigui- 
ties on the leading-term perturbative expansion of order A 2 /m 2 from the ultraviolet 
renormalon sequence associated to the divergence of the perturbative series for the 
spectral function P~5|-fn|]. The present note is mainly devoted to a quantitative dis- 
cussion of the impact of UV renormalon ambiguities on the determination of a s (m 2 ). 
On the basis of the accumulated knowledge on renormalon behaviour, we address the 
question of what is the theoretical error on a s (m 2 ) and examine possible ways to 
decrease it. We discuss the magnitude of UV renormalon corrections obtained from 



explicit calculations [p!q -|2Tl], which although based on unrealistic simplified schemes, 
are nevertheless interesting to some extent. We then study a number of improved 
approximants for the perturbative series, based on a change of variable in the Borel 
representation |i4| , |T7|| , such as to displace the leading UV renormalon singularity at a 
larger distance from the origin than the first infrared (IR) renormalon. The spread of 
the resulting values of a s (m 2 ) obtained by different approximants, at different renor- 
malisation scales is exhibited as a measure of the underlying ambiguities. Finally, on 
the basis of mathematical models, we discuss the prospects of an actual improvement, 
given the signs and the magnitudes of the computed coefficients, the size of a s (m 2 ) 
and what is known of the asymptotic properties of the series. Our conclusion is that 
a realistic extimate of the theoretical error cannot go below 5a s (m 2 ) ~ ±0.060, or 
5a s (m 2 z ) ~ ±0.006. 

The organisation of this article is as follows. In sect. 2 we summarise the basic 
formulae and discuss different procedures to do the integration over the circle that 
differ by resumming or not an infinite series of "large 7r 2 terms". We discuss the 
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relative merits of the various procedures and their scale dependence. In sect. 3 we 
introduce the problems related to the divergence of the perturbative series, we review 
the Borel transform method and the renormalon singularities. In sect. 4 we derive 
some useful formulae obtained in the Borel space after integration on the circle. In 
sect. 5 we consider the explicit form for the leading UV renormalon singularity derived 
in perturbation theory in the large Nf limit, Nf being the number of flavours. This 
limit is not meant to be realistic, but, for orientation, we evaluate the quantitative 
impact that such an UV renormalon would have on the determination of a s (m*). 
We find that this effect is rather small. In sects. 6, 7, which contain the main original 
results of this work, we introduce and study a number of improved approximants that 
could in principle suppress the ambiguity from the leading UV renormalon. We study 
the combined effects of different, a priori equivalent, procedures, different accelerators 
of convergence and different choices of the renormalisation scale. We also study in 
a simple mathematical model under which conditions for the known coefficients of 
the series the accelerator method leads to a better approximation of the true result. 
Finally, in sect. 8 we present our conclusion. 



2. Basic Formulae and Truncation Ambiguities 

The quantity of interest is the integral over the hadronic squared mass s in r 
decay of a function R(s) analogous to R e + e -(s), weighted by a phase-space factor. In 
the limit of massless u, d, s quarks we have 



J o m A T \ m A T J \ J 
R{s) is proportional to the imaginary part of a current-current correlator: 

N N 

R(s) = — Im n(s) = — : [n(s + ie) - U(s - ie)} . (2.2) 

The normalization factor N is defined in such a way that, in zeroth order in pertur- 
bation theory, R(s) = 3. In turn, the correlator Tl(s) is related to the Adler function 
D T (s), defined in such a way as to remove a constant: 

DJs) = -StMI(s). (2.3) 
ds 
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By first integrating by parts and then using the Cauchy theorem one obtains for R T 
the result 



2m J\a\=m% s 



ds ( s 



mi 



1 + DM 



mi 



(2.4) 



The Adler function D T has a perturbative expansion of the form: 

oo 

D T (s) = D° T Y, D n a(-s) n ~ D° [l + D ia (-s) + D 2 a 2 (-s) + D 3 a 3 {-s) + 



n=0 



(2.5) 

where a = a s /ir, D® = 3(1 + 5) where 6 is a known small electroweak correction, and, 
for Nf = 3 in the MS scheme, 



D a = l 
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(2.6) 



where ((3) = 1.20206 and C(5) = 1.03693, C A = N c = 3 and C F = (JV£-l)/(2jV< 



c 



4/3. The quantity (3 = (11CU — 2iV/)/12 is the first beta function coefficient 

/5(a(/i 2 )); /3(a) = -/3a 2 (l + /3'a + . . .) (2.7) 



,da(/i 2 ) 



d/1 2 



(/3 = 9/4, /?' = 16/9 for iV/ = 3). 

The expansion in eq. ( |2.5|) defines the Adler function at all complex s with a cut 
for s > 0. In the spacelike region, where s < 0, a(— s) is real and given asymptotically 
by Gu 2 > 0): 



1 



-+ y + /J log=£=/ 3 log=£. (2. f 



a( — s) alW) /r 
If we want o(— s) at some complex value of the argument, e.g. s — —\s\ exp(i#), we 
can use the formula 

a(-s) = (2.9) 



l + pa(\s\)i9 
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where the angle 6 is —n on the upper tip of the cut for s real and positive, +ir on the 
lower tip and zero on the negative real axis. The more accurate two-loop expression 
is given by 

a(|s|) 

a(_S) = 1 + (3a{\s\)id + P'a(\s\) log(l + (3a{\s\)id) ' (2 ' 10) 

The expansion for R(s) (for s real and positive) can be obtained from that of D T 
by the relation 

^-kivjjw- (2 - n) 

Performing the integration by using the expansion for D T in eq. 
sion in eq. (|2.10 ) for a complex argument, one obtains 



5|) and the expres- 



and 



R(s) = D° T (l + Fio(s) + F 2 a 2 (s) + F 3 a 3 (s) + ... 



Fx = Di; 



D 2 ; 



D, 



[2.12) 



(2.13) 



The origin of the /3 2 7r 2 /3 term is easily understood. By using the one-loop expansion 
for a(s), eq. ( p.9|) , one gets 



R{s) = Dl 



1 + 



1 



log 



1 + in(3a(s) 



+ ... 



p 2 7T 2 a 3 ' 



1 + a(s) 



+ ... 



2n(3i ° 1 — m(3a(s) 

(2.14) 

We have the following observations on this result. First, for Nf = 3, the coefficients 
Fi 2,3 in the expansion of R(s) for r decay coincide with those of R e + e - because the 
potentially different terms proportional to (J2Qi) 2 , with Qi being the quark charges, 
vanish in this case. Second, we observe that eq. Q2.13D is obtained by a truncation 
of higher-order terms in the quantity f5 2 a 2 7T 2 ~ 0.7 with a = a(m 2 ). Note that a 
similar problem of truncation arises when the integration over the circle in eq. (|2.4j ) 
is performed. In the early treatments of this problem (e.g. in ref. [§]) the expression 
of a(—s), which appears on the circle, is taken from eq. ( p.lOj ) and expanded in a 
consistently to the order a 3 . With this procedure, one obtains 

^(bnp) = D o ^ + ffia ( m 2) + H2a ^ m l) + H 3 a 3 (m 2 T ) + ...), (2.15) 



with 



H x = 1; H 2 = 5.2023; H 3 = 26.3666. 



(2.16) 
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We will refer to this result as BNP formula (for the authors of ref. ||). More recently 
in ref. @| it was advocated that a better procedure for performing the integration 
on the circle is to keep the full three- loop expression for a( — \s\e te ), according to the 
formula 

i?( LP ) = — (l-z) 3 (l+z) 1 + D 1 a(-zm 2 T ) + D 2 a\-zm 2 T ) + D 3 a\-zm 2 T ) + . 



2ni J\z\=i z 

(2.17) 

It is this procedure which is currently adopted (LeDiberder-Pich, or LP method). At 
fixed experimental value of R T , the two procedures lead to values of a s (m 2 ) that differ 
by terms of order 5a s {m 2 ) ~ ([3 2 a 2 ir 2 )a 2 ~ 0.01 or ((3 2 a 2 n 2 ) 2 a ~ 0.05. These "large 
7r 2 terms" always arise when one goes from the spacelike to the timelike region (e.g. 



similar terms arise [fj^] when one relates Drell-Yan processes to electroproduction or 
fragmentation functions to structure functions). There have been many discussions 
in the past on the opportunity of resumming these terms ||23|| . If there was a good 



argument to consider the expansion for D T (s) in some respect superior to that for 
R(s) it could be worthwhile at low energies to keep the expression in eq. ( |2.14| ) in 



its resummed form rather than to expand in ft 2 a 2 TT 2 . A glance at eqs. ( |2.6| ) shows 
that, for n < 3, there are no explicit "large 7r 2 terms" in the coefficients D n of the 
expansion for D T . Similarly when the integration on the circle is performed with the 
complete formula for a(— \s\e l6 ) all terms are kept up to order a 3 (f3 2 a 2 ir 2 ) n , i.e. up 
to order a 3 one expands in a but keeps /3 2 a 2 n 2 unexpanded. Sure enough this sequel 
of terms exists in reality, so why not take them into account? However, the counter- 
argument is that there are in perturbation theory terms involving ir 2 that arise from 
origins other than the spacelike-timelike connection and, in any case, there are many 
terms of the same general magnitude (for example, the term proportional to /3 2 in 
D 3 , eq. (|2]6|)), so that the advantage of keeping this particular class of terms is likely 
to be completely illusory. 

One can further consider the scale dependence of the different procedures. One 
can write 

jj>(BNP) = D o ^ + &i fl(/i2) + H 2 a 2 ^ 2 ) + # 3 aV) +■■■), (2-18) 

where 
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H 3 ( f i 2 ) = H 3 -2H 2 (3\og^ + H l (3 2 log 2 -f -0/3' log . (2.19) 




Analogously, we can study the scale dependence in the LP method. In this case we 
have 



where the D n coefficients are related to the D n as the H n to the H n in eqs. ( |2.19| ). 

The results of the LP and BNP methods are shown in fig. where, assuming a 
measured value of 3.6 for R T , we show the corresponding determination of a s (m 2 ) 
as a function of the renormalization scale /x. For comparison, we also show (dashed 
curve) the determination obtained with a simplified LP method, in which we use the 
one-loop expression eq. ( |2.9|) for a. For /x > 1 GeV there is less /x dependence if the 
resummed expressions are used. This stability is often taken as a possible indication 
that resumming is better. 

In conclusion, we agree that the resummed formulae provide a less ambiguous 
result with respect to a change of scale than the unresummed expression. However, 
it is true that a priori it is not possible to guarantee that a more accurate result is 
obtained in one way or the other. As a consequence the spread shown in fig. [I] for 
different choices of /x and of procedure is to be taken as a real ambiguity. In particular 
the large discrepancy at /x ~ 1 GeV is a genuine signal of trouble, especially in view 
of the fact that several proposed scale-fixing procedures lead to small values of /x, e.g. 
minimal sensitivity, BLM scheme etc. ||24j| . A small value of \x is also suggested by 
physical considerations, because the average hadronic mass is well below m T . So, on 
the one hand, one cannot sensibly reject the option of small values of /x. On the other 
hand, the corresponding value of aim 2 .) becomes very ambiguous at small /x. 

3. Renormalons and Borel Transformation 




2iTi J\z\=i z 




3 (l+z) [l + D^i-z/j, 2 ) + D 2 a 2 (-zfi 2 ) + D 3 a 3 (-zfi 2 ) + ...], 

(2.20) 



In the following discussion we first study the properties of D T (s) in itself, and only 
later we consider the integration over the circle. The problem that we now consider 
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Figure 1: The value of a s {m^) obtained from R T = 3.6, as a function of the 
renormalization scale fi. 
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is the well known fact that the series for D T is divergent. Indeed one can identify 
sequences of diagrams, depicted in fig. called "renormalon" terms [ 13| - pl[| , that 
provide the leading behaviour at large n for the n-th coefficient of the expansion for 
D T . The renormalon contribution is of the form: 




+ 2 



n bubbles 



Figure 2: Renormalon diagrams for a two-point correlator. 



D n ~C k nW k (^j [l + 0(l/n)] (n large). (3.1) 

Note the n! behaviour which implies that the series is divergent. Here, C k and 7^ > 
are not known for the real theory but only, to some extent, in the large- Nf expansion 
18|-|TJ , and the index k runs over a discrete set of values: 



k = —1, —2, —3, . . . Ultraviolet (UV) Renormalons 

k = +1(?), +2, +3, . . . Infrared (IR) Renormalons (3.2) 

In the above list we omit the contribution from instantons, which appear at rather 



large values of k. They have been computed in ref. and shown to be small. The 
UV or IR renormalons arise from the limits of large or small virtuality, respectively, for 
the exchanged gluon(s). As hinted by the question mark, the k = +1 IR renormalon 



is probably absent in perturbation theory, but the issue is not really settled |T6| - pT| 1. 
The absence of this term is necessary for the consistency of the assumption that all 
non-perturbative effects can be absorbed in the condensates. In the following we will 
assume that the k — 1 IR, renormalon is indeed absent. 
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In view of the divergence of the perturbative expansion, one can possibly give a 
meaning to the quantity 

d(a) = ^-l = D ia + D 2 a 2 + D 3 a 3 + ... (3.3) 



by the Borel transform method J27 ] . One defines the perturbative expansion of the 



Borel transform B(b) of d(a) by removing the n! factors: 

oo n l2 

B(b) = £ D n+1 — = Di + D 2 b + D 3 — + . . . (3.4) 

Then, formally 

POO 

d(a) = / dbe- b ' a B{b) (3.5) 



o 

in the sense that the expansion for B(b) reproduces the expansion for d(a) term by 
term. What is needed for d(a) to be well defined is that the integral converges (this 
cannot be true at all s [|T~3| , [T~7[| because of the singularities of d(a) in the s plane, but 
this problem can be neglected in our context) and that B(b) has no singularities in 
the integration range. But, as already mentioned, the large-n expansion of B(b) leads 
to singularities on the real axis. In fact, at large n, B(b) is essentially given by a 
geometric series: 

B(b) ~ C k J2^ k (j^J ~ C k Y( lk + + less singular terms (3.6) 

so that it is singular at b — k//3. Thus the UV renormalons correspond to singularities 

at b = —1/(3, —2/(3, ■ ■ ■ and the IR renormalons at +2//?, +3//3, As a consequence, 

the convergence radius of the expansion for B(b) near the 6-origin is determined by 
the UV renormalon at b = — 1//3, independent of the existence of the IR renormalon 
at b = +1/(3. 

Thus, the perturbative expansion for B(b) can be directly used only to perform 
the integration up to b = +1/(3. The contribution from b = +1/(3 up to b = oo, 
where the expansion is not valid, could typically lead to terms of order A 2 /s (or even 
worse). For example, if B(b) is sufficiently well behaved at b = +1/(3 and at b = oo, 

/•oo 

Ad(a) = / dbe- b/a B(b) ~ aB(l//3)exp(-l/0a) ~ aB(l/(3)A 2 /s, (3.7) 
Ji/p 
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where we used a -1 ~ /31og(s/A 2 ) and the fact that the exponential cuts away all 
large-6 contributions so that B(b) was approximated by its value near b = +1/(3. 
From a different point of view, at large n, the series for d(a) is dominated by the UV 
renormalon behaviour with D n ~ rCnl^— (3) n . At fixed small a, the individual terms 
\D n \a n first decrease with n, then flatten out and eventually increase because of the 
n\ factor. The best estimate of the sum is obtained by stopping at the minimum, 
for n ~ n opt , given by \D n \a n ~ |D„_i|a n_1 , or n opt ~ l//3a. From the theory of 
asymptotic series [p^j , the corresponding uncertainty Ad(a) is of order \D nopt \a n ° pt : 

\D n Ja^ ~ {l/(3af{l/(3a)\{(3af' Pa) 

~ (l//3a) 7 (l//?a) (1//3a) e WMy/2n/l3a(l3a)WM 

~ (l//3a) 7 e (1//3a) v /27r//3a ~ A 2 /s x logarithms, (3.8) 



(where the Stirling approximation was used: n! ~ n n e~ n y/27rn). Thus one could 
improve the accuracy of the perturbative expansion by computing more subleading 
terms until n ~ n opt is reached and then add a residual term of order A 2 /s. Note that 
the estimate n opt ~ l//3a ~ 4 indicates a rather small value. However this estimate is 
obtained from the behaviour of the leading UV renormalon series, while there is no 
alternation of signs and in general no evidence of renormalon behaviour in the few 
known terms of the series. 

While an accuracy of order A 2 /s is what one gets in practice from the three- 
loop expression of d(a), it is true that, in principle, if there is no IR renormalon 
at b = +1//3, d(a) can be better defined. In fact, as the location of the leading 
UV renormalon at b = —1/(3 is not in the integration range, there is the possibility 
of defining B(b) by analytic continuation up to b = +2/(3. If this is realised then 
the remaining ambiguity, of order (A 2 /s) 2 , is unavoidable because the corresponding 
singularity at b = +2/(3 is on the real axis, so that an arbitrary procedure to go around 
it must be defined and the difference between two such procedures would be of that 
order. However, since operators of dimension 4 do exist in the operator expansion, 



this ambiguity can be reabsorbed in the non-perturbative condensate terms ||lT3l4T5f| . 
We also understand that the absence of the IR renormalon at b = +1//3 is necessary 
for the consistency of the SVZ approach because operators of dimension 2 are 
absent in this channel. If there would be a singularity at b — +1/(3 the corresponding 
ambiguity could not be absorbed in a condensate. There are indeed indications in 
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perturbation theory that the first IR renormalon does not appear [T^j2T|. However 
there could be non-perturbative sources of breaking of the operator expansion at non- 
leading level. After all no theory of confinement could be built up from perturbation 
theory and renormalons. But, in practice, independent of the existence of the IR 
renormalon at b = +1/(3, the accuracy to be expected from the first three terms in 
the expansion for d(a), as they have been used so far in the actual determination of 



a 



is of order A 2 /s. 



4. The Integration over the Circle in the Borel Transform 
Formalism 



In this section we show that the integration over the circle in eq. (|2.4j ) for R T is 
particularly simple in the Borel representation. Starting from eqs. ( |2.9| ) and ( |3.5| ) we 
have 

We work in the approximation where the two-loop coefficient f3' in the beta function 
is neglected. Then, according to eq. 02.90 , we can replace l/a(—s) by l/a + ij39, where 
a is a{\s\) = a(mj), invert the integration order and write s/m 2 = — exp(i#): 

r= / dbe~ b/a B(b)— tde(l + e ie ) 3 {l-e ie )e- ib(3e . (4.2) 
Jo 2ni J-tt 

The integration is easily performed, with the result 

Jo K) [3b{l3b-l){(3b-'Z){(3b-A)<K 

= / dbe~ b/a B{b)F{f3b). (4.3) 
Jo 

We see that, in first approximation, the effect of going from a(—s) to a(\s\) by inte- 
grating over the circle is to multiply the Borel transform B(b) by the factor F(/3b). 
For real x, the function F(x) is shown in fig. |3|. It is an entire function in the whole 
complex plane with a good behaviour at infinity on the real axis. The factor F(x) has 
simple zeros at the location of all UV renormalons and also of all IR renormalons with 
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Figure 3: Plot of the function F(x). 
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the exception of those at (3b = 1 (if any), 3 and 4. Since, in general, the correspond- 
ing singularities are not simple poles, they are not eliminated, but their strength is 
attenuated (this point will be discussed in more detail in sect. 5.). Equation (|4.3| ) 
obviously coincides with the LP approach 0] in the limit (3' = (3" = 0. One can also 
repeat the procedure by expanding in a(/i 2 ) instead of a(mj), according to eq. fl2.20|) . 
The resulting fi dependence is shown in fig. [I] together with the analogous results for 
the BNP and the LP formulae. 



5. Large- Nf Evaluation of Renormalons and their Resum- 
mation 

As well known, the typical renormalon diagrams of QED and QCD can be eval- 
uated in the large- Nf limit and their structure is simple in this limit |T8| - pT| . In the 
abelian case the large- Nf limit corresponds to the large-/? limit, and this is believed to 
be true also in the non-abelian gauge theory in spite of the fact that the beta function 
in this case cannot be evaluated only in terms of vacuum polarisation diagrams. The 
sequence of dominant terms in (3 n generalises the terms in (3 and (3 2 that appear in 
D2 and D3 respectively. As argued in a recent paper, ref. [Q, the determination 



of the exact behaviour of the UV renormalon series may be very different from the 
one indicated in the large-/? limit. While we do not know the complete form of the 
leading UV renormalon at b = — 1/(3, we can nevertheless compute, for orientation, 
the quantitative impact of its approximate form at large (3 on the determination of 
a(mj). From eq. (42) of ref. ]2T|] one obtains the large-/? expression of the contribution 
of the leading UV renormalon at b = —1/(3 to the Borel transform Bib): 



B{b) = V 5 / 3 £ [7 + 2n] (-x) n = V 5 / 3 



7 lr 



1 + x 



1 + x' 



(5.1) 



where x = (3b and the factor e~ 5 ^ 3 transforms the result from the MOM into the 
MS scheme. We observe that the leading UV renormalon is a double pole. Since 
in the large-/? limit (3' can be neglected, the corresponding expression for R T in the 
approximation of eq. ( |4.3| ) is appropriate, and it turns the double pole into a simple 
pole. 

We now study the numerical effect of including the whole renormalon series with 
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respect to a truncated result up to the order b 2 . We first consider the impact on 
d(a), i.e. before the integration on the circle. In order to get the correction to d(a) 
from the higher-order terms in the UV renormalon we must subtract from B(b) its 
expansion up to 0(b 2 ) and perform the inverse Borel transform, eq. (|3.5|) : 



Ad{a) 



1 

P Jo 



dx e 



-x//3a 



B(x) 



e - 5/3 (7-9x + llx 2 



9 



(5.2) 



For (3 = 27/12 = 2.25 and a = 0.12 (or pa = 0.27) one finds Ad 3.7 x 10~ 3 

which corresponds to a ~2% increase in the value of a s (m 2 ) at fixed d(0.12) = 0.155 
(5a s (m 2 T ) ~ 0.007). 

We now repeat the same exercise for the function r, given in eq. (O), obtained 
after integration over the circle. We compute the variation 



Ar(a) = - 



roo r 9 

/ dxe~ x/Pa B(x) - -e" 5/3 (7-9a; + llx 2 
jo L 9 



FfaO 



(5.3) 



Numerically we find Ar(0.12) ~ -6.0 x 10" 3 which, at fixed r(0.12) = 0.220, 
again corresponds to 5a s {m 2 ) ~ 0.007. Thus the extra factor F(pb) has practically 
no influence on the effect on a s (m 2 ) of the nearest UV singularity. 

In conclusion the overall effect of the UV renormalon singularity in this model is 
small and not much changed by the integration over the circle. 



6. Search for More Convergent Approximants 

Assuming that indeed there is no IR renormalon at b = +1/P one can in principle 
try to obtain by analytic continuation a definition of the Borel transform, valid on the 
positive real b axis up to b = +2/P, outside the radius of convergence of its expansion. 
We now discuss how the analytic continuation could be implemented in practice. 

Starting from eq. (|3.4j ) we can make a change of variable |14| , |17| z = z(b) with 
inverse b = b(z), z(0) = and z(oo) = 1 (so that the interval from to oo in b is 
mapped into the to 1 range in z), such that the IR singularities are mapped onto 
the interval between zq = z(2/P) and 1 and the UV singularities are pushed away at 
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\z\ > zq. Changing variable one obtains 





db 


/ dz 


lo 


dz 



e' b(z)/a B(b(z)). 



Using the expansion 



the series 



goes into 



b(z) = ciz + c 2 z 2 + . . . 



Bib) = D ± + D 2 b + D 3 - + ... 



B{b{z)) =D ± + D 2 c lZ + {D 2 c 2 + D 3 ^)z 2 + 



(6.1) 

(6.2) 
(6.3) 
(6.4) 



which is convergent up to z = Zq, while the original b expansion was convergent only 
up to b = corresponding to z(l/(3) < z . The improved approximation for d(a) 
is therefore given by 



rzo 

d(a) ~ / 

J z 



dz 



db 


e -b{z)/a 


dz 





r 2//3 , i 

/ dbe- b/a 
Jo 



D 1 + D 2Cl z + [D 2 c 2 + D 3 ^)z 2 + 



D l + D 2Cl z(b) + [D 2 c 2 + D 3 ^)z(b) 2 + 



(6.5) 



where the full expression of z as function of b is inserted in the integral. In this way, 
an infinite sequence of terms is added to the b expansion. For a small, the upper limit 
of integration can be replaced with infinity without significant effect. 

One possible example is given by [17]: 



z(b) 



^TTpb+1 



b(z) 



Az 

/3(1-*) 2 ' 



(6.6) 



In this case the first UV singularity is at z — — 1, and all higher UV renormalons are 
on the unit circle \z\ — 1. IR renormalons are between zq = (v^3 — l)/(\/3 + 1) and 
z = 1. In this example, c\ — 4//3, c 2 = 8//3. Other examples are 



z(b) 



(3b 



k + (5b 



b{z) 



kz 



(3(1 -zY 



(6.7) 



with k — 1,2 or 3. Also in these cases the first IR renormalon at b — 2/(3 becomes 
the closest singularity to z — 0, while the UV are pushed further away. Here we have 
ci = c 2 = k/(3. 
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Before discussing numerical applications, we observe that the present method 
relies simply on the position of the IR and UV renormalon singularities in the Borel 
plane and not on the nature and the strength of the singularities. We have seen 
that the integration over the circle in eq. (|4.3|) does not change the position of the 
singularities in the b plane, but simply affects their strength. Thus, we can as well 
consider the effect of the accelerators on the expansions for i?(. BNP ) given in eq. ( 2.18|) 
or on the LP expression of eq. ( [2.20|) . For example, the improved version of eq. 
simply becomes 



/ dbe~ b/a 
Jo 



c 2 

D x + D 2 c lZ (b) + (D 2 c 2 + D 3 ^)z(b) 2 + 



F{(3b). (6. 



We now consider the following numerical exercise. We assume that experiments 
have measured R T = 3.6. We then compute a(m 2 ) with the LP formula as a function 
of the scale /i, and we perform the same calculation applying our acceleration proce- 
dures to the LP method The results are shown in fig. |j. We see that relatively 
large differences in the fitted value of a s {m 2 ) are obtained, especially at large \x for 
different accelerators and in comparison to the non-accelerated formulae. We do not 
see a priori compelling reasons to prefer one or the other procedure. The fact that a 
priori equivalent methods lead to results with a sizeable spread must be considered 
as an indication of a real ambiguity. Even if we only consider the method of ref. [|J 
for the integration over the circle, it is impossible to go below an uncertainty of the 
order Sa s (m 2 ) ~ ±0.050 for /i in the range from 1 to 3 GeV. The ambiguity becomes 
even larger if we extend the comparison to the formulae with truncation in f3 2 a 2 n 2 
a-la BNP (fig. §). 



7. Study of More Convergent Approximants in a Model 

The method for accelerating the convergence discussed in the previous section 
only relies on the position of the singularities of the Borel transform and not on their 
nature and strength. It looks rather surprising that one can compensate for the effect 
of renormalons without actually knowing their form in detail. In this section we study 
a simple mathematical model to clarify under which conditions the method can be 
successful, in the sense that it provides a better approximation to the true result. 
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Figure 4: Effect of the accelerators on the determination of a s (m^.) with the 
LP method, for R T = 3.6. The curves b, c and d refer to the change of variable 
of eq. ( [6, ID , while e refers to eq. ( 6.t ). 
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We consider as a model the case where the Borel transform is exactly specified by 



b 2 00 / n + 1 — 1 « 
B true (b) = l + D 2 b + D 3 ^ + pJ2 [ n \ i-(3b) n 



n=3 



1 + D 2 b + D 3 - + p 



(1 + ^6)t 



- 1 + 7/36 - 



7(7 + V ) 2 



• (7-1) 



The added sum stands for the higher-order contribution that could arise from a leading 
U V renormalon at b = — 1 / (3 with a degree of singularity specified by 7 and a fixed 
overall strength given by p. We will take p = 1 in the following discussion. It is 
convenient to re-express eq. ( pM|) in terms of 2 = (3b: 



B true (x) = 1 + D 2 x + D 3 — + 



1 



(1 + ar)' 



7(7 + 1) 2 
1 + 72 x 



(7.2) 



where -D 2 = -D2//3 and D 3 = D 3 /(3 2 . Similarly we can introduce B pert (x) and 
B aC cei{ x ), the perturbative Borel functions without and with acceleration, respec- 
tively: 

(7.3) 



x 



B p ert{x) = 1 + D 2 X + D 3 — 



t-2 



Baccel(x) = 1 + D 2 C lZ (x) + (D 2 C 2 + D 3 ^-) Z (x) 2 , 



(7.4) 



where Ci j2 = (3cx )2 . In all cases the corresponding d function, d true , d pert and d acce [ is 
given by 

Ma) = / dxe~ x/l3a B(x). (7.5) 
Jo 



We consider the ratio 

dtrue d acce i 



H 



dtrue dp er t 



l-D 2 I 2 {(3a)-D 3 I 3 {(3a), 



(7.6) 



where the quantities J 2 ,3 are given by 

HP a ) = TTo~\ i dxe~ x/Pa (c 1 z(x) + c 2 z{ 
lo((3a) Jo v 

i . v 

dxe~ x/(3a (cjz{x) 2 - x 2 ) 



x) — X 



I 3 ((3a) 



1 1 



2 I Q ((3a) Jo 



(7.7) 
(7.8) 



-21- 



and 

roc 

Uf3a) = / dxe- x/l3t 
Jo 

Clearly, \H\ < 1 is the condition for the acceleration method to be successful. In 
particular for H = d trU e and d acce i coincide. For each value of (3a and 7, in a given 
model specified by z{x) and the corresponding coefficients Ci 2 , the condition H = 
is satisfied on a straight line in the plane D2,D%, while the inequality \H\ < 1 is 
satisfied in a band defined by two straight lines parallel to the H = line. In figs.|6|-|8] 
we plot the lines H = for (3a = 0.27 for fixed 7 and z(x) given by eq. (|6.6| ) (case 
labelled by 0), or by eq. ( |6.7| ) with k = 1 or 2 (cases 1 and 2). The values of 7 in 

2.0 
1.5 
1.0 

|p 0.5 
0.0 
-0.5 

-2 -1.5 -1 -0.5 0.5 1 

D 2 

Figure 6: Lines corresponding to H = /or i/ie method "0" (solid), "1" 
(dashed) and "2" (dotted) for 7 = 0.5. The circle corresponds to expansion of 

1/(1 +X)T. 

figs. |6|-[8] are 7 = 0.5, 1, 2. We see that for each choice of z(x) the lines have different 
negative slopes. The lines tend to cross each other in a region of the plane not far 
from the point D 2 = —7, -D3 = 7(7 + 1), i.e. the values that correspond to the first 
few terms of the expansion of the asymptotic function 1/(1 + x) 1 . The region where 



[1 + x) 



1 + 7X 



7(7+1) 2 

x 



(7.9) 





Figure 7: .4s in Fig. |6|, for 7 = 1. 
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Figure 8: As in Fig. |6|, for 7 = 2. 



-24- 



the lines cross is more sharply defined if 7 is small, i.e. if the asymptotic function is 
not too singular. In figs. ^| p7o| we show the bands \H\ < 1 for 7 = 0.5, 1 in cases 
and 1. 



4 



3 



2 







-1 



1 1 1 1 1 1 1 
I 7=0.5 \ \ 


111 

\ ^^^^ 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


\ ^ - 

\ 

\ - 

1 ^ 
1 1 1 1 1 1 1 1 



-4 -3-2-10 1 2 



D 2 

Figure 9: Bands corresponding to \H\ < 1 for the method "0" (solid), "1 
(dashed) for 7 = 0.5. 



The conclusion is that the method for accelerating the convergence works well only 
if the coefficients D 2 , D 3 resemble those of the asymptotic series, in other words if the 
known terms in the expansion are sufficiently representative of the asymptotic series. 
In particular we see that it is very unlikely to get an improvement if the coefficients 
D2 and D3 are of the same sign, as is unfortunately the case for the series of interest 
for us (see eq.6). 

From a different point of view we now consider the simple function B(f3b) given 

by 

B{f3b) = - l — (7.10) 
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Figure 10: As in Fig. [9|, for 7 = 1. 
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and we plot the relation of the exact result 



POO 

Ma) = / dxe~ xlf * a B[ 
Jo 



X 



(7.11) 



with its accelerated or non-accelerated series approximants, as a function of (3a and 
of the order of the expansion. The results obtained for the accelerating function z(b) 
given in eq. (|6.6|) (the case labeled by 0) and 7 = 1,2 are shown in figs. ITT ,|12. We 




Figure 11: Effect of the resummation technique described in the text, for a 
function with Borel transform B(b) = 1/(1 + /36) 7 ,7 = 1 



see that when, as in this case, the coefficients of the expansion coincide with their 
asymptotic form the accelerated formulae provide a much better approximation to 
the true result, more so if the singularity is weaker (i.e. 7 is smaller). The non- 
accelerated formulae are only good for small enough /3a. The physically interesting 
case of f3a ~ 0.27 appears to be at the limit of the range where the non-accelerated 
formulae are acceptable. 




Figure 12: Effect of the resummation technique described in the text, for a 
function with Borel transform B(b) = 1/(1 + /36) 7 ,7 = 2 
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8. Conclusion 



The determination of a s {m T ) or a s (m z ) from r decay is nominally very precise 
and the experimental errors are extremely small at LEP. Certainly the dominant 
ambiguity is at present the theoretical error. The nominal precision is large because 
in the massless limit no explicit 1/m, 2 corrective terms are present in the operator 
expansion. But it has become clear by now that one cannot sensibly talk of power- 



suppressed corrections if the ambiguities in the leading term are not under control [25 



We think that there is no real theorem that prevents non-perturbative corrections in 
the coefficient function of the leading term in the operator expansion at the level of 
A 2 / m 2 . Equivalently, there could be an IR renormalon singularity at b — which 
would create an irreducible (being located on the integration path) ambiguity of order 
A 2 /m 2 . In all-order perturbative evaluations of the singularity pattern in the Borel 
plane the IR singularity at b = +1//3 is probably absent, a result consistent with the 
idea that all irreducible ambiguities can be reabsorbed in condensates. Even if the 
IR renormalon singularity at b = +1//3 is indeed absent, the radius of convergence of 
the expansion is limited by the leading UV renormalon singularity at b = — 1/(3. If 
this disease is not cured or cannot be cured the resulting ambiguity is still of order 
A 2 /m 2 . In principle the problem could be solved if the exact nature and strength of 
the singularity was known, by simply taking its effect into account in the evaluation 
of a s (m 2 ), along the way indicated in section 5 in the case of the estimate of the 
singularity in the unrealistic limit of large (3. But the exact determination of the 
singularity appears to be beyond the scope of presently known methods. In the 
actual case where the UV renormalon singularity at b — — 1/(3 is not specified, one 
can still try, in principle, to bypass the problem by a transformation of variables 
that pushes the leading UV singularity to a larger distance from the origin than 
the first IR renormalon singularity at b = +2/(3. Expanding in the new variable is 
equivalent to add a specified infinite sequel of terms to the original expansion. The 
convergence of the series should be improved by these accelerators of convergence and 
the ambiguity decreased. We have studied the quantitative effect of implementing a 
number of such accelerators with different choices of the renormalisation scale /i. An 
indication of the size of the ambiguities on a s (m 2 ) is obtained from the spread of the 
results for different starting formulae (e.g. with a(—s) taken in the integration over 
the circle in its renormalisation group improved form or in a fixed order truncated 
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expansion), different accelerator functions and different choices of the renormalisation 
scale. The difference between the resummed or truncated expression for a(—s) on the 
circle are expecially large at small values of //, while the variations induced by the 
different accelerators are expecially pronounced at large values of fi. The relative 
stability of the unaccelareted result of ref. || versus changes of \i appears as largely 
accidental in that the accelerated formulae based on it are much less stable at large 
(jl. It was argued in ref. ]19) that if one expands in a s (/i 2 ) instead of expanding 
in a s (m 2 ) the scale dependence of the UV renormalon correction becomes of order 
(A 2 /m^)(mJ//i 2 ) 2 . Can then one be safe if /i is chosen sufficiently large? Clearly 
in the true result the sum of the perturbative terms plus the remainder must be 
scale independent. When \i is changed, the number of terms to be added before 
the series becomes asymptotic changes and must compensate for the difference. The 
increased sensitivity of the accelerated formulae at large /i is not encouraging for 
invoking that large [i is safer. All together, from fig. 4 we find it difficult to imagine 
that the theoretical error on the strong coupling can be taken smaller then, say, 
8a s {m1) ~ 0.050 (which approximately corresponds to 5a s (m 2 z ) ~ 0.005). 

The accelerator method is based on the mere knowledge of the position of the 
singularity and not on its precise form. Clearly such a method can only work if the 
known terms of the expansion carry enough information on the asymptotic form of 
the series. We have quantitatively confirmed this statement by studying the problem 
on a simple mathematical model where the true result is known. The performance of 
different accelerators is studied as a function of the coefficients of the first few terms. 
These results indicate that there is little hope of improving the ambiguity from the 
leading UV renormalon because the first few coefficients of the actual expansion show 
no evidence for the asymptotic behaviour, in particular no sign alternance. This last 
argument (as well as the one on the /x dependence of the UV renormalon) can be 
interpreted in different ways. If one is a great optimist, he can argue that the series 
does not resemble at all to the renormalon asymptotics, hence the normalisation 
of the renormalon term is very small (as is the case for the explicit form of the 
singularity obtained in the large Nj limit). Or, if one is more cautious, as one should 
be in estimating errors, he can say that since the known terms do not show sign 
of asymptotia, they are dominated by subasymptotic effects and cannot be used to 
estimate the remainder. In this spirit we do not propose the accelerators as a better 
way to determine the true result but simply as a criterium to evaluate the theoretical 
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error. In fact, while all accelerators tend to increase the resulting value of a s , the 
amount of the upward shift is sizeably different for different accelerators. 

In order to bypass all possible objections one should be able to fit at the same time 
a s (m 2 ) and C2, the coefficient of 1/m 2 . corrections to R T . Note that in the ALEPH 
moment analysis || C2 is fixed to zero while the coefficients of some higher-dimension 
operators are fitted. This is not very relevant to the main issue. If C2 is not fixed it is 
found that the sensitivity to a s (m 2 ) is much reduced. In an interesting paper Narison 
p9| attempted to put an upper bound on C2 from the data on e + e~ — > hadrons. This 
is an important issue that would deserve further study. Our interpretation of the 
analysis of ref. |29[ is that values of C2 of order (500 MeV) 2 are not at all excluded. 
Narison |29| derives a more stringent limit \C 2 \ < (374 MeV) 2 but we feel he relies 
too much on the so called optimisation procedure. Indeed, something that should 
be a constant in a dummy variable turns out to be a steep parabola. The value at 
the tip is taken, with a small error, as the best estimate because of the vanishing of 
the derivative at that point, instead of considering the span of the results in a priori 



reasonable range for the irrelevant parameter. In a recent paper |30| an estimate of 
C2 from Argus data on hadronic r decay was obtained and the results are compatible 
with |C 2 | < (500 MeV) 2 . 

We ignored here other possible sources of error beyond those arising from higher 
orders in perturbation theory. These include errors from the freezing mechanism for 
a s , errors from the translation of a s (m 2 ) in terms of a s (m|), from the region of the 
circle integration near the positive real axis and so on. These errors are presumably 
smaller || than our current estimate of the error from higher order terms in the 
perturbative espansion. Taking all the other uncertainties into account we end up 
with a total theoretical error around 8a s {m 2 z ) ~ 0.006. As a result, in spite of the 
fact that our estimate of the error is larger than usually quoted, the determination of 
a s (m 2 z ) from r remains one of the best determinations of the strong coupling constant. 
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